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ABSTRACT 

The effects of the UV background radiation on the formation of sub-galactic clouds are 
studied by means of one-dimensional hydrodynamical simulations. The radiative transfer of 
the ionizing photons due to the absorption by HI, Hel and Hell, neglecting the emission, is 
explicitly taken into account. We find that the complete suppression of collapse occurs for the 
clouds with circular velocities typically in the range Vc ^ 15 — 40 km s~^ and the 50% reduction 
in the cooled gas mass with ~ 20 — 55 km s^^. These values depend most sensitively on 
the collapse epoch of the cloud, the shape of the UV spectrum, and the evolution of the UV 
intensity. Compared to the optically thin case, previously investigated by Thoul & Weinberg 
(1996), the absorption of the external UV photon by the intervening medium systematically 
lowers the above threshold values by AVc ~ 5 km s^^. Whether the gas can contract or 
keeps expanding is roughly determined by the balance between the gravitational force and the 
thermal pressure gradient when it is maximally exposed to the external UV flux. Based on 
our simulation results, we discuss a number of implications on galaxy formation, cosmic star 
formation history, and the observations of quasar absorption lines. In Appendix, we derive 
analytical formulae for the photoionization coefhcients and heating rates, which incorporate the 
frequency /direction-dependent transfer of external photons. 
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1. Introduction 

Photoionization of primordial gas is well known to have great impacts on galaxy formation and the 
thermal history of the universe. The observations of the QSO absorption spectra (Gunn & Peterson 1965) 
imply that the intergalactic gas has been highly ionized before redshift z ^ 5. In fact, the existence of 
an intense ultraviolet (UV) background radiation that can photoionize the universe is inferred from the 
proximity effect of Lya forest ai z — 2 ^ A (e.g. Baljtlik, Duncan & Ostriker 1988; Bechtold 1994; Giallongo 
et al. 1996). The origin of the UV background radiation could be attributed to the emission from QSO's 
and/or young galaxies formed at high redshifts (e.g. Couchman 1985; Miralda-Escude & Ostriker 1990; 
Fukugita & Kawasaki 1994). 

The UV background radiation, once produced, largely affects the subsequent formation of structures. 
In particular, the formation of low-mass objects is suppressed via photoionization and heating associated 
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with it (Umcmura & Ikcuchi 1985; Ikcuchi 1986; Rccs 1986; Couchman & Rees 1986; Ikeuchi, Murakami 
& Rcos 1988, 1989; Efstathiou 1992; Babul & Rccs 1992; Chiba & Nath 1994; Babul & Ferguson 1996; 
Okoshi & Ikeuchi 1996; Haiman, Rees, & Loeb 1997; Kepner, Babul & Spergel 1997, among others). Such 
suppression may partly resolve some of the shortcomings of the hierarchical models of structure formation, 
such as an excess number of faint galaxies predicted than actually observed (e.g. White & Frcnk 1991; 
Kauffmann, White & Guiderdoni 1993; Cole et al. 1994). In view of this, several authors have simulated 
the thermal and dynamical evolution of the intergalactic medium under the UV background (Umemura & 
Ikeuchi 1984; Bond, Szalay & Silk 1988; Murakami & Ikeuchi 1990, 1993; Cen & Ostriker 1992; Zhang, 
Anninos & Norman 1995; Quinn, Katz & Efstathiou 1996; Thoul & Weinberg 1996; Weinberg, Hernquist & 
Katz 1997; Navarro & Steinmetz 1997). For instance, Thoul & Weinberg (1996) concluded that the objects 
with circular velocities 1^ ^ 30 km s~^ are prevented from collapsing at 2 ^ 2: ^ 5. These previous studies, 
however, were based upon the assumption that the medium is optically thin against the ionizing photons. 
Gas clouds in reality become optically thick in the course of contraction and it is by no means trivial how 
much the cloud evolution is altered when the radiative transfer of ionizing photons is incorporated. In 
addition, only the limited range of collapse redshift has been explored previously. As the mean density 
of the universe depends strongly on redshift in proportion to (1 + z)^, the effects of the UV background 
radiation should also vary, with an increasing importance of the radiative transfer at higher redshifts. The 
possible evolution of the UV background intensity may have further impacts on the cloud dynamics at 
different redshifts. 

In this paper, we study the effects of the UV background radiation on the dynamical evolution 
of primordial clouds, by means of onc-dimensional, spherically symmetric hydro dynamical simulations, 
incorporating the radiative transfer of the ionizing photons. We use analytical approximations of the 
photoionization coefficients and heating rates derived in Appendix, which take explicit account of the 
frequency /direction-dependent radiative transfer due to the absorption by HI, Hel and Hell. They are 
simply expressed as a hmction of column densities of each species from the cloiid boundary and are 
applicable to the medium with an arbitrary density profile. By using these analytical formulae, we have 
achieved a significant reduction of the computational time and are thus able to explore broad ranges in 
cloud mass scale, collapse redshift, and parameters of the external UV field. The results are compared 
quantitatively with those of optically thin calculations. Based on the results of numerical simulations, we 
further predict the global production rate of cooled gas in the standard Cold Dark Matter (CDM) universe 
and discuss its implications on galaxy formation and cosmic star formation history. We also discuss the 
imprints of the UV background on the quasar absorption lines, such as the observability of the helium Lya 
forests, which can be tested by future observations. 

The plan of this paper is as follows. In §2, we describe the numerical model used in this paper. §3 
presents the results of our simulation and §4 is devoted to the discussion. Finally §4 summarizes our 
conclusions. Wherever necessary, the following cosmological parameters are assumed for definiteness; the 
density parameter flo = 1, the Hubble constant h = i?o/(100km Mpc""'^) = 0.5, the baryon density 
parameter flh — 0.1, and the amplitude of the density fluctuations erg = 0.6. The value of (Tg is chosen to 
match the observed local abundance of galaxy clusters (e.g. Viana & Liddle 1996; Eke, Cole & Frenk 1996; 
Kitayama & Suto 1997). 



2. Method 
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2.1. Basic Equations 

We simulate the dynamics of a spherically sjonmetric bound system exposed to the diffuse UV 
background radiation. The system is a mixture of baryonic gas and coUisionless dark matter, with the mass 
ratio oi fib : Qq — Qb = 1 : 9. The evolution of these components is described by the following equations: 
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Here r, m, p, P, T, u, and are the radius, mass, density, pressure, temperature, internal energy per unit 
mass, and mean molecular weight in units of the proton mass nip, respectively. M(< r) is the total mass 
interior to r, 7i and C arc the heating and cooling rates per unit volume, G is the gravitational constant, 
and fcfi is the Boltzmann constant. Wherever necessary, the subscripts b and d denote baryon and dark 
matter, respectively. The adiabatic index is fixed at 7 = 5/3 throughout the paper. 

The above equations are solved using the second-order-accurate Lagrangian finite-difference scheme 

described in Bowers & Wilson (1991) and Thoul & Weinberg (1995). The shocks arc treated with the 
artificial viscosity technique (Richtmycr & Morton 1967; Umcmura 1993). The shells are binned equally in 
mass and their numbers are Nb = 500 for baryonic gas and = 5000 for dark matter. We have performed 
runs with {Nb,Nci) = (2000,20000) and (300,3000) and confirmed that our results are robust against the 
changes in the resolution. We have also checked that our code reproduces accurately the similarity solutions 
for the adiabatic accretion of coUisional gas and for the pressure-less collapse onto a black hole (Bertschinger 
1985). The total energy of the system is conserved by better than a few percent in any runs reported in this 
paper. 



2.2. Radiative Processes 

The baryonic gas is assumed to have the primordial composition with hydrogen and helium mass 
fractions X = 0.76 and Y = 0.24, respectively. At each time-step, starting out from the cloud boundary 
into the interior, we successively solve for the ionization equilibrium among photoionization, collisional 
ionization and recombination, together with the penetration of the external UV field, as described in 
detail later in this section. We then compute the heating/cooling rates due to photoionization, collisional 
ionization/excitation, recombination, thermal bremsstrahlung, and Compton scattering with the cosmic 
microwave background radiation. Unless otherwise stated, we use the rates and coefficients summarized in 
Fukugita & Kawasaki (1994), which corrects a few typos in Cen (1992). 
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The external UV field is taken to be isotropic and have the power-law spectrum: 

Jin(j/, z) = J2i(z) ( — ^ X lO^^^erg s-^ cni-^ str"! Hz"!, (6) 

where J2i(z) is the intensity (in proper coordinates) at the Lyman limit of hydrogen {hi/m = 13.6eV) and 
a is the spectral index. Observations of the proximity effect in the Lya forest suggest J21 — 10*°'^ at 
z = 1.7 - 4.1 (Baljtlik, Duncan & Ostriker 1988; Giallongo et al. 1996; Cooke, Espey & Carswell 1997; 
Savaglio et al. 1997), but its value is still highly uncertain at other redshifts. Theoretical predictions of 
the value of J21 needed to reionize the universe at high redshifts range from J21 ~ 0.1 to even 100 (e.g. 
Fukugita & Kawasaki 1994; Gnedin & Ostriker 1997; Haiman & Loeb 1998; Madau, Haardt & Rees 1999). 
The predicted epoch at which reionization occurs also has a large uncertainty between z ~ 50 and 6. In the 
present paper, unless otherwise stated explicitly, we fix the onset of the UV background at zyv = 20, and 
consider the following four cases at z < zuv- 

1. J21 — 1 and a = 1, 

2. J21 — 0.1 and a = 1, 

3. J21 — 1 and a = 5, 

4. Evolving J21 and a = I, where 

r (i+£)-6 6<z<zw 
J21 = < 1 3 < z < 6 (7) 

[ i^)' 0<z<3. 



The values a ^ 1 and 5 are chosen to mimic the spectra of quasars and massive stars, respectively. By the 
above form of the UV evolution (eq. |^ ) , we attempt to study the effects of a late reionization and a decline 
of the UV intensity at low redshift. 

The incident UV spectrum is modified due to the radiative transfer as it penetrates into the gas cloud. 
In the present paper, we explicitly take account of the frequency /direction-dependent absorption by HI, 



Hel, and Hell, using equations (A15) and (AI6) in Appendix. Having performed the frequency and angular 
integrations of the transfer equation neglecting the emission term, the analytical formulae are obtained for 
the photoionization coefficients and heating rates in the plane-parallel slab. They are simply expressed as 
a function of column densities of individual species measured from the boundary, and lead to a significant 
reduction in the computational time. By means of these formulae, we solve simultaneously for the ionization 
equilibrium in each gas shell and the UV radiation field processed between the outer boundary and that 
shell. The ionization state of the gas and the radiation field are computed iteratively until the abundances 
of HI, Hel and Hell all converge to the precision better than 1%. 

Note that the above treatment of the radiative transfer systematically overestimates the effect of 
absorption in the following respects; 1) since the solution of the transfer equation in the plane-parallel 
geometry is applied to a spherical cloud, the photon path lengths from the outer boundary are overestimated, 

2) at each radial point, ionizing photons coming from the inner 27r steradian of the cloud is ignored, and 

3) the emission (or scatter) of incident photons is neglected. It should thus provide a conservative limit in 
which an external UV field is maximally weakened and has the minimal effects on cloud evolution. This 
is in fact complementary to a conventional approximation of the optically thin medium (e.g. Thoul & 
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Weinberg 1996), in which the external UV field is likely to have the maximal effects on the cloud dynamics. 
In what follows, therefore, we examine both of these limiting cases and perform quantitative comparisons 
between them. In so doing, we hope to bracket the true answer, which is still very difficult to solve in a 
fully self-consistent manner. 

The approximation of the ionization equilibrium is correct if the recombination time-scale ircc is 
shorter than the dynamical time-scale t^yn = l/\/Gp. Assuming a high degree of ionization, this yields the 
following condition for the electron density Ue (e.g. Vedel, Hellsten & Sommer-Larsen 1994): 

> 7.1 X lO-^Ti V-^V^) cm-3, (8) 



1 + X J \ 10 

where T4 = (T/10^ K), and we have adopted troc = (aH«e)~^ using the recombination rate for hydrogen 



an given in equation (A27). The above condition is satisfied in most situations considered in this paper 



except at the outer envelopes of clouds at low redshifts. For example, in a cloud collapsing at z = 3, 



typically the out-most ^ 5% of the gas has the density below the above value at z = 3 (see §2.3 for the 
definition of the collapse epoch). For these shells, the assumption of the ionization equilibrium will result 
in overestimating the fractions of neutral species and hence in overestimating the effect of absorption even 
further. In practice, however, the physical abundances of neutral species in these shells are negligibly small. 
Thus the approximation of the ionization equilibrium is well justified in the present analysis. 



2.3. Initial Conditions 

We start the simulations when the overdensity of a cloud is still in the linear regime. The initial 
overdensity profile is assumed to have the form of a single spherical Fourier mode for both baryon and dark 
matter components: 

c / N c /^N sin(fcr) 

S^ir)^S,{0)^-^, (9) 

where k is the comoving wavenumber, and 6^(0) is the central overdensity. Throughout the paper, we fix 
Si(0) = 0.2. Assuming that the initial perturbation is dominated by the growing mode, the initial velocity 
profile is given by 

v,{r)^H,r(l-^-^^), (10) 



where Hi is the Hubble parameter at the initial epoch, and 6{< r) is the volume averaged overdensity 
within radius r. The outer boundary is taken at the first minimum of Si{r), i.e. kr = 4.4934, at which 
^(< r) vanishes and the shell initially expands at the speed Hir. As in Haiman, Thoul & Loeb (1996), we 
define the baryonic mass enclosed within this radius as the bound mass Af bound, and that enclosed within 
the first zero of 5i{r), i.e. kr — tt, as the cloud mass Mdoud- These masses are related to each other via 
0.342[l-H 0.304(5i(0)]Mb ouncl 

= 0.363Mb ound- 

Having fixed the initial density profile, we vary the initial redshift Zi and the cloud mass Mdoud to 
simulate different collapse epochs and mass scales, i.e., circular velocities or virial temperatures. We define 
the central collapse redshift Zco and the cloud collapse redshift Zc respectively as the epochs at which the 
inner-most gas shell and the shell enclosing Afcioud would collapse to the center in the absence of thermal 
pressure. The circular velocity Vc and the virial temperature T^ir of the cloud are defined using A/doud and 
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clS 

J..,, . ,^o9xio'(^)(^||=;^)^",i + ,, K (12) 

Unless otherwise stated explicitly, we study three combinations of Zi, Zco, and listed in Table ^. In 
this table, we also list the turn-around redshifts corresponding to Zco and Zc, ztao = 2^/'^(l + Zco) ~ 1 and 
zta = 2^/'^(l + Zc) — 1, respectively. In the low-redshift collapse, the gas is exposed to the UV background 
radiation from the linear regime. In the middle-redshift collapse, the cloud center is close to turn-around 
at the onset of the UV radiation zyv = 20. In the high-redshift collapse, the cloud center has already 
collapsed prior to zuv- 

At an initial epoch, the gas is assumed to have the uniform temperature given by 

/ 2.726(1 + 2.) ifz.>200 

I 548(^)^'^-^^ ifz,< 200, ^^'^ 

taking account of the fact that the matter is tightly coupled to the cosmic microwave background at z ^ 200 
(e.g. Anninos & Norman 1994; Chieze, Teyssier & Alimi 1997). 



2.4. Boundary Conditions 

If the gas is able to lose sufficient energy by radiative cooling, it falls towards the center nearly at the 
free-fall rate. As the density increases near the center, the cooling time-scale ^cooi = wp/ \H — L\ and the 
dynamical time-scale ^dyn become extremely small, and infinitely large number of time-steps are required 
in the simulation. In order to avoid this, we introduce the following criteria. If a gas shell reaches below 
some minimum radius rmin caid satisfy tcooi/^dyn < 0.01, then we regard it as having cooled and collapsed; 
we move the shell to the center and ignore in the rest of the simulation, except in the calculation of the 
gravitational force. When dropping a shell to the center, we assume that the next shell expands inward 
adiabatically to cover the volume of the dropped shell. The evolution of each shell is traced until both of 
the above criteria are fulfilled. We adopt as rmi„ the radius at which the system attains rotational support 
(e.g. Padmanabham 1993): 

where rta is the turn-around radius of the gas shell, and Ata is the dimensionless spin parameter attained 
by turn-around. We fix Ata = 0.05 based on the results of the numerical simulations (Efstathiou & Jones 
1979; Barns & Efstathiou 1987; Warren et al. 1992). Equation ( [T^ ) also gives the scale under which 
the assumption of spherical symmetry breaks down. The baryon mass that has collapsed to the center is 
denoted by A/coii- 

The dynamical time-scale also becomes very small when collisionless dark matter shells approach the 
center. We therefore rebounce the dark matter shells near the center (Spitzer & Hart 1971; Gott 1975). 
As noted by Thoul & Weinberg (1995), it is desirable to choose the rebouncing radius rj-cb that is small 
enough not to alter the evolution of gas shells, but not too small to degrade energy conservation. In order 
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to achieve these, we set r^cb equal to of a gas shell that encloses 0.01 — O.lMbound depending on a run. 
Throughout each run, Ti-ob is fixed at this value. 

At the outer boundary, we adopt the free boundary condition, i.e. the pressure outside the cloud Pout 
is zero. We have also performed runs with a mirror boundary condition (Pout = P/Vb, where Pn^ is the 
pressure of the outermost gas shell) and confirmed that our results remain essentially invariant. 



3. Numerical Results 

3.1. Dynamical evolution 

Figure |l] illustrates the trajectories of gas shell radii in the low-redshift collapse (zc — 0.5) for a cloud 
with Vc — 32 km s^^ (Mdoud — 9 x 10^ Mq). In the absence of the UV background, the gas shells initially 
comoving with the Hubble expansion will eventually turn around due to gravity and start to contract. As 
the virial temperature of the cloud is Tvir ~ 4 x 10^ K, the gas can lose energy efficiently by radiative cooling 
and collapse towards the center. If the UV background radiation is present, on the contrary, the evolution 
of the cloud is significantly altered. Figure |^ (a) shows that the cloud is prohibited to collapse completely 
under a constant UV with J21 = 1 and a = 1. Once, however, the UV amplitude is reduced to J21 = 0.1 
and the absorption is taken into account, the inner part of the cloud may start to collapse (Fig. |l|b). 
Alternatively, a softer spectrum a — 5 enables larger number of shells to collapse by z = (Fig. |l|c). The 
cloud collapse is also promoted if the UV background radiation evolves according to equation ^ (Fig- |l|d). 

In order to study the above mentioned features in more detail, we plot in Figure |^ the evolution of the 
temperature T, HI fraction X^ii = »t.hi/('^hi + J^hii), HI column density measured from the outer boundary 
A^Hi, and Hell fraction AThcII = ''t-hoIi/("-HcI + '^HoII + 't-HoIIi) in the runs shown in Figure |l|. For J21 = 1 and 
a — 1, temperature rises to T 10^ K and all the gas is photoionized almost instantaneously at the onset 
of the UV background zyv = 20 (Fig. ||a). At this epoch, the cloud density is still almost equal to the mean 
of the universe (nn ~ 2 x 10~^ cm~^) and the external UV field can penetrate into the cloud center, even if 
the absorption is considered. The thermal pressure of the heated gas takes over the gravity and the cloud 
keeps expanding. 

For a smaller UV amplitude J21 =0.1, Figure ||(b) shows that the cloud evolves very differently 
with/without the absorption. If the absorption is considered, the central part of the cloud is not ionized 
instantaneously, but stays neutral until the gas density decreases significantly due to expansion. It is 
noticeable that the temperature still rises almost instantaneously to T ~ 10"* K. This is because the 
time-scale for heating is much shorter than that for ionization (e.g. Gnedin & Ostriker 1997). One can 
roughly estimate the ratio of these time-scales in the case of a = 1 as 

^hcat _ 3kBT /p-l 

iion 2'Hhi / ™ 

_ j 0.285r4 optically thin limit, 

[ 7.28 X IQ-^T^T^I^^ optically thick limit, 

where thi — 6.3 x 10~^®(A^Hi/cm^^) is the HI optical depth at the Lyman limit z/hi, and Fhi and 
TYhi are photoionization coefficient and heating rate derived for the pure hydrogen gas in Appendix 
(eqs [1A22| -[ |A25| ). At zuv, the HI column density at the center is A^hi ^ 10^° cm"^ in Figure |(b), which 



gives the optical depth of thi = 630 and thcat/^ion = 8.5 x 10 in the above equation. For a softer 
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spectrum, this ratio becomes larger and the heating becomes relatively slower as illustrated in Figure ||(c). 
In this case, the pressure gradient between the cold neutral center of the cloud and the hot ionized envelope 
causes to push the gas slightly inward near the cloud center at z ~ 10 (Fig. |^ c). 

Figure ^(b) further indicates that as the gravity overcomes the pressure support, the gas shells can turn 
around and the density starts to rise again. The gas is first heated by shock and adiabatic compression, 
but it will start to lose energy and shrink to the center once the radiative cooling becomes efficient. 
Consequently the neutral core is formed at the center. The absorption of the external UV field promotes 
the cloud collapse and the formation of the neutral core, by both reducing the number of photoionized ions 
and lowering the heating rate. 

The cloud evolution is very sensitive to the spectral index of the UV background. For a softer spectrum 
a — 5, there are smaller numbers of high energy photons, and hence larger numbers of Hell, than for a ^ 1 
(Fig. 1^). The larger number of shells are thus able to collapse within a given time for a — 5. 

If the UV background evolves according to equation (0), the gas in the cloud is ionized gradually as in 
Figure ||(d). As the UV intensity drops at 2; < 3, the neutral fraction starts to rise again and the inner 
shells start to collapse to the center. 

Figure || shows the radial profiles at z = in some of the runs discussed above. The shock front (if 
exists) is at r ~ 10^ kpc. The density profile inside the shock front agrees well with the self-similar solution 
p (X r~^-^^ (Bertschinger 1985). The central part of the cloud is self-shielded against the external UV 
photons and the gas becomes neutral. In the case of (J2i,a) = (0.1,1) with absorption, the simulated 
profiles of Xhi and A^hi roughly agree with an analytical estimation for the pure hydrogen gas presented 



in Appendix. For instance, inserting into equation (A29) the values read off from the figure at r ~ 2 



kpc, (nn ~ 10~^cm~^, T ^ lO'^K, and Am ~ 0.1), we obtain thi ^ 2, i.e. A^m ~ 3 x 10^''cm~^, in good 
accordance with Figure ^c). On the contrary, such an estimation breaks down for the softer spectrum 
(a — 5) because of the larger abundances of Hel and Hell. 

The effects of the UV background radiation become less significant as the cloud mass increases. Figure 
^(a) shows that the gas collapse is still delayed by the UV radiation for a cloud with Vc = 50 km s^^ 
(Afcioud — 3 X 10^ M0, Tvir ~ 9 X 10^ K). The evolution, however, is only slightly altered for a larger cloud 
with Vc = 100 km s'^ (Mdoud ^ 3 x 10^° Mg, Tvir ~ 4 x 10^ K). The effects of absorption, compared to the 
optically thin case, also become insignificant at this mass scale as the collisional ionization dominates the 
photoionization. 

Figures H and || illustrate the trajectories of gas shells in the middle-redshift collapse {zc — 3). Since 
the density of a cloud is higher than the low-redshift collapses discussed above, the larger amount of gas is 
able to collapse for a given circular velocity Vc = 32 km s~^. At the onset of the UV background zyv = 20, 
the cloud center is about to turn around, and the density is njj ^ lO^'^ cm^'^, roughly a factor of 5 larger 
than the mean density of the Universe. At this density, the inner part of the cloud can be kept self-shielded 
against ionizing photons until it collapses if the absorption is considered (Fig. ^,c,d). Even in such cases, 
the gas is still heated promptly to T ~ lO'^K for a = 1, because of the high efficiency of heating process 
over ionization mentioned above (eq. fl^]). 

The effects of the UV background radiation is further weakened if one goes to even higher redshifts. 
Figure |^ shows the gas shell trajectories in the high-redshift collapse {zc — 10). At zuv, about 10% of 
-^^bound have already collapsed, and a number of shells either have turned around or are about to turn 
around. In addition, Compton cooling further promotes the gas collapse at z ^ 7. 
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3.2. Critical mass-scales for the collapse 



In order to quantify the results presented in §3.1, we plot in Figure H the fraction of gas mass which 



has cooled and collapsed to the center after the UV onset zyy = 20 as a function of circular velocity. 
The output epochs correspond to t = OMc, tc, 2tc, where tc = io/(l + z^Y^"^ and to is the present age 
of the universe. In the low and middle redshift collapses {zc = 0.5 and 3), there exists a sharp cut-off in 
the collapsed fraction and this threshold lies in the range Vc = 20 ~ 50 km s~^ depending on different 
assumptions on the UV background. In the high-rcdshift collapse {zc — 10), on the contrary, the threshold 
is not so distinct and smaller cloud is able to collapse against the external UV. 

We define the cut-off velocity U^f"* and the half-reduction velocity V^^^'^, respectively, as the velocities 
at which the mass collapsed between zuv and Zc is 0.05 and 0.5 of that with no UV. Figure ^ illustrates 
that these quantities depend most largely on the UV spectral index a, the evolution of the UV flux, and 
the collapsing redshift Zc (we have added to this figure the results for the runs with Zc = 0, 1.4, 4.8 and 
6.2). Compared to the optically thin cases, the adoption of the radiative transfer systematically lowers V^^'^ 
and y^^^'^ by about 5 km s~^. 

For a constant UV flux (Figs ^ a-c), both V^^'^ and V^^^^ decrease with increasing z^, because of higher 
gas density and stronger Compton cooling. At Zc ^ 7, the central part of the cloud begins to collapse 
prior to zyv and V^^^ falls below 15 km s~^, corresponding to Tvir ~ 8000K (note that l/*^"' and V^^'^^^ are 
evaluated from the amount of gas cooled after the UV onset). We simply omit the data points from the 
figures at such low velocities, since it is beyond the scope of our current framework to study these regimes 
where the molecular cooling dominates. 



In the presence of the UV evolution (Fig. gd), Vj,"^"' and V"^^^'^ rise rapidly with time at high redshifts 
and then begin to drop at low redshifts. However, the peaks of these quantities are at Zc ^ 1.5, which is 
much later than the peak of the UV flux (z = 3 — 6). This is because the kinetic energy of expanding gas 
particles delays the cloud collapse signiflcantly even after the UV flux starts to decline. The late infall of 
gas after the decline of the UV background is therefore suppressed until Zc ^ I. 

So far we have flxed the onset of the UV background at zuv ~ 20. This choice is rather arbitrary given 
the large uncertainties in the reionization history of the universe. A possibility of much later reionization 
(smaller zuv) can be covered effectively by the UV evolution of the form (^. In order to examine an 
alternative possibility of much earlier reionization, we performed a few runs with zuv = 50 and added the 
resultant 1/^"' to Figure ||(a). The differences from the fiducial zuv — 20 runs are significant only for the 
high-redshift collapse {zc ~ 10), because the gas shells infalling at z = 20 — 50 are mainly affected. At such 
high redshifts, the effect of the radiative transfer of UV photons becomes large because the cloud density is 
very high. 

Figure ^ shows the baryon mass M^^^^ corresponding to Vj?^'^ as a function of Zc. Also plotted for 
reference are the Jeans mass with T = lO^K and the masses corresponding to 1,2,3a density perturbations 
in the standard CDM universe with Hq — 1, h — 0.5, ilb = 0.1, and erg = 0.6. The cloud well above 
Tvir = lO'^K is significantly affected by the UV background radiation. For a constant UV flux, M^^^^^ rises 
at a rate faster than (1 -I- Zc)~^/^, which is expected for the Jeans mass with constant temperature. This is 
because Vj^"*'^ increases with decreasing Zc as shown in Figures ^ (a)-(c). For the evolving UV, M^^^^ keeps 
growing even after the decline of the UV background and the growth slows down considerably only at z ^ 1. 

Figure |l^ further indicates that the collapse of ^ 2a density fluctuations in the standard CDM 
universe is highly suppressed at low masses by the UV background radiation. For example, the cooled mass 
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from a lo" fluctuation is reduced to less than 50% at Mdoud ^ 10^ Mq and Zc /t 2 for (J21, a) = (1, 1), and 
at il/cioud 2 X 10^ Mq and Zc }t 3 for {J21, a) = (1, 5) (see also §4.2 for more discussion on the implications 
on structure formation). 



4. Discussion 

4.1. Comparison with analytical estimates 

The results of our simulations are compared with analytical estimates on the density-temperature 
diagrams (e.g. Efstathiou 1992; Chiba & Nath 1994) in Figure |l^. We divide the simulation results into 
three categories and denote their virial density ng"^ = 3.80 x 10~^(ri6/i^/0.025)(X/0.76)(l + Zc)^ cm~^ and 
the virial temperature Tvir by different symbols; 1) efficient cooling (circles), if the gas mass Mdoud can 
cool before tc + ^dyn, 2) moderate cooling (triangles), if Mdoud can cool only within the present age of the 
universe io, and 3) inefficient cooling (crosses), if otherwise. For comparison, three analytical relations are 
computed, using the same UV parameters as the simulations except for assuming optically thin in all cases; 
icooi = ^dyn (solid line), tcooi = to — tc (dashed), and Tvir = 7oq (dotted), where Teq is the temperature at 
which the cooling and heating rates balance {Ti = C). The former two analytical relations are evaluated 
using at Zc, while the last one using = '^h'/S at the turn- around epoch zta- 

Figure ^ shows that our simulation results under a constant UV are in reasonable agreement with 
the analytic estimates defined above; 1) efficient cooling lies in the region where tcooi < ^dyn, ^cooi < to — tc 
and Tvir > ?eq, 2) moderate cooling mainly lies where idyn < ^cooi < to — tc and Tvir > Toq, 3) inefficient 
cooling lies where icooi > to — tc or Tvir < Teq. Such correspondences degrade near the boundaries of each 
region. Compared to the optically thin runs, cooling efficiency is systematically enhanced to some extent if 
the absorption is taken into account. 

It should be noted that, in order to achieve above agreements, it is essential to evaluate the relation 
Tvir = Teq at Zta, not at Zc (if evaluated at Zc, the dotted lines in Fig. |ll| are shifted downwards by a factor 
of 8). This is because the cloud evolution largely depends on the photoionization prior to the collapse. 
Roughly speaking, whether the gas can contract or keeps expanding is determined by the balance between 
the gravitational force and the thermal pressure gradient when it is "maximally exposed" to the external 
UV flux, i.e., when the gas attains the maximum value of J2i/7iH- For the gas exposed to the constant UV 
flux from the linear regime, this corresponds to the turn around. The above statement also applies in the 
presence of the evolution in the UV intensity. 



4.2. Cosmological implications 

The suppression of low-mass objects by the UV background radiation has profound implications on 
cosmology and galaxy formation. Given the fraction of cooled gas in objects of different masses and collapse 
redshifts from the simulations, we can deduce the abundance distribution of collapsing objects in the entire 
Universe. Specifically, we compute the baryon mass density that cools and collapses per unit redshift as a 
function of the collapsed mass Mcoii of baryons in a halo and the collapse epoch Zc as 

dMd-z = dMd-z ' 
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where /coii is the coUapse fraction of a cloud relative to the pressure-less case, and (PN^^i^{M,z)/dMdz 
is the comoving number density of halos with mass M that collapse and form per unit redshift at z. We 
use the latter quantity rather than a conventional mass function, e.g. that of Press & Schechter (1974), 
because f^oii obtained in our simulations is more closely related to clouds just collapsing at a given redshift. 
Unfortunately, there is not yet a fully self-consistent analytical formalism to compute this quantity, while 
several approaches have been proposed (e.g. Bond et al. 1991; Lacey & Cole 1993; Sasaki 1994; Kitayama 
& Suto 1996a, b; Manrique & Salvador-Sole 1996). In this paper, as a working hypothesis, we adopt the 
halo formation rate given by equation (15) of Kitayama & Suto (1996a) with the threshold mass M{ ~ M/2 
(see their paper for detail). We have checked that our results are qualitatively unchanged by adopting an 
alternative approach by Sasaki (1994). 

Figure |l^ illustrates the baryon mass distribution of collapsing objects in the standard CDM universe. 
In the absence of pressure, collapsing objects have a broad mass distribution with an increasing fraction 
of high mass objects at lower redshifts. Once the gas pressure is taken into account, high mass objects 
are reduced by the low cooling efficiency, and low mass objects by the UV background. As a result, the 
abundance distribution at 0c ~ 3 agrees well with the observed mass range of galaxies. Note that the low 
mass end of the distribution has a tail, not a sharp cut-off. This is because the objects just above V]^"' can 
have a wide range of A/coii for almost identical Afcioud- Thus the collapse of dwarf-sized objects of mass 
Mco\i ^ 10^ Mq is still possible at Zc ^ 0.5, but such objects are expected to have lower baryon fraction 
than normal galaxies. This point can be tested by future observations. 

One can integrate the above distribution over mass to obtain the total amount of baryon that is 



collapsing at a given epoch. Figure 13 shows the evolution of this quantity (denoted by "cold gas -t- star"). 
We have also deduced similarly from our simulations the baryon mass density that becomes bound by the 
gravity of clouds but not yet cooled at Zc- Such a component, denoted by "hot gas", is mainly contained in 
objects corresponding to groups and clusters of galaxies. We further define "unbound gas" as the gas that 
becomes unbound due to the UV background but would have collapsed at Zc in the absence of pressure. 



Figure 13 indicates that the production rate (per Hubble time at a given epoch) of "cold gas + star" has a 
peak at Zc ~ 2 — 3 and it is rather insensitive to the UV parameters. This seems to suggest a high efficiency 
of star formation activity at these redshifts. On the other hand, the production rate of "hot gas" and that 
of "unbound gas" simply increases and decreases with time, respectively (see also Barkana & Loeb 1999; 
Cen & Ostriker 1999). 

Figure |lj further compares the predicted production rate (per year) of "cold gas + star" in the 
standard CDM universe with the observed cosmic star formation rate (SFR, Madau et al. 1996, 1998; 
Lilly et al. 1996; Connolly et al. 1997) compiled by Totani et al. (1997) and Totani (1999, private 
communication). This comparison is only demonstrative, since the predicted curves simply set an upper 
limit to the global SFR in the universe and their normalization can be shifted in proportion to the baryon 
density parameter Nonetheless, it is noticeable that the predicted production rates of "cold gas -I- 
star" under the UV background show a steep rise at z ^ I, in good qualitative accordance with what is 
actually seen in the SFR data. At higher redshifts, they continue to rise up to z ~ 4 and then flatten. This 
is somewhat different from a rather sharp decline at z ~ 2 — 4 in the observed SFR. It should be kept in 
mind, however, that these SFR data are still very uncertain due to the dust extinction and may be shifted 
upward by large factors (e.g. Pcttini et al. 1998). Recent detections of star forming galaxies in the submm 
band also suggest rather high SFR of ^ 2 x 10~^/i Moyr^^Mpc"^ at z 2 - 4 (Small, Ivison & Blain 
1997; Hughes et al 1998; Berger et al. 1998). More detailed discussion on this subject should await higher 
precision data from future infrared observations. 
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4.3. Imprints on quasar absorption lines 

The effects of the UV background radiation on the intergalactic medium can be observed most directly 
by means of absorption lines of quasar spectra. In order to discuss the impacts of current results on such 
observations, we plot in Figure |l^ the line of sight column densities through the cloud Nj {j =HI, Hel, 
Hell) as a function of impact parameter p. Here, we restrict the observability of these elements to the 
column densities > lO^'^cm"^, and define the critical impact parameter pcrt within which the column 
density is greater than lO^'^cm^^. As the column density of each element is too small to be detected 
except at the central region in the case of J21 = 1 and a = 1, we focus on the other two cases of UV 



parameters shown in Figure 15, In Tables]^ and |3|, we summarize Pcrt and the ratio of its square value, 
which corresponds to the ratio of area, i.e., the relative expected number. These tables suggest that a large 
number of helium forest lines are expected especially under the soft UV spectrum. Recent observations by 
HST in fact indicate the detections of numerous Hell lines (Hogan, Anderson & Rugers 1997; Reimers et 
al.l997; Anderson et al. 1999). 

If we fit the lines in Figure |l^ at Nj > lO^'^cm^^ by a single power-law Nj cx p^", then the column 
density distribution of each absorption line is written as dAf /dNj cx (Nj)^^ with /3 = (n + 2)/n. In Table 
^, we summarize n and (3 values. Since these are crude numbers, we merely present the values in the case 
of evolving J21 at Zc = 0.5 and ignore the dependence on redshift and UV parameters. At iV^j > lO^^cm"^, 
the column density distribution of Hel is similar to that of HI while that of Hell is a little steeper. These 
points can be checked by future observations. 

The above results are also clearly seen in Figure |l^, which plots the line of sight column densities 
of Hel and Hell against the HI column density. The relations among these column densities are rather 
insensitive to Vc and Zc, and are summarized in Table |^ (values are given for = 32 km and Zc = 0.5). 
In order for the Hel lines to be detected at N^^^ > lO^'^cm"^, the corresponding HI forest lines are required 
to have TV^j > 10^* — lO^^cm"^. On the other hand. Hell forest lines are more easily detectable and are 
even saturated at N^^ > lO^^cm^^. 

In addition to the high observability of the helium absorption lines at the UV wavelengths, our results 
further suggest an interesting feature in the HI absorption lines at high column densities. When the 
line of sight passes near the central region of a cloud that has just collapsed, the gas is expected to be 
neutral but the temperature can be as high as lO^K, because the time-scale of ionization is longer than 
that of the photoionization heating (eq. ||l5|). It should be kept in mind that if the hydrogen molecules, 
ignored in the present paper, are present, they could allow cooling below 10'*K. In any case, the existence 
of the warm neutral gas could be checked by the HI absorption lines with a large Doppler parameter and 
the 21cm emission line with high spin temperature. Incidentally the recent observations of damped Lya 
systems (Lane et al. 1998; Chengalur & Kanekar 1999) suggest the detection of neutral gas with the spin 
temperature ~ lO^K, which is much higher than what is found in normal spiral galaxies. 



5. Conclusions 

We have shown that the formation of the sub-galactic clouds is greatly prevented by the UV 
background radiation even if the transfer of the external UV photon is taken into account. Within the 
range of parameters investigated in this paper (J21 < 1, a = 1, 5, Zc = — 10), the complete suppression of 
collapse occurs for the clouds with circular velocities typically in the range Vc '--^ 15 — 40 km s^^ and the 
50% reduction in the cooled gas mass with Vc ^ 20 — 55 km s~^. These values depend most sensitively on 
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the collapse redshift Zc and the slope of the UV spectrum a. 

The evolution of the UV background also affects the above thresholds in a significant manner. The 
decline of the UV intensity at z <, 3 can decrease the threshold circular velocities at lower redshifts. This 
effect, however, is delayed until z ^ 1.5 due to the kinetic energy of gas particles attained at higher redshifts. 
In fact, whether the gas can contract or keeps expanding is roughly determined by the balance between the 
gravitational force and the thermal pressure gradient when the gas attains the maximum value of J2i/nu- 

As far as the evolution of the gas down to T ~ lO'^K is considered, the radiative transfer of the ionizing 
photons has the moderate effect. Compared to the optically thin case, the absorption of the UV photons 
by the intervening medium systematically lowers the above threshold values by AVc ~ 5 km s~^. Once the 
evolution below T ~ lO^K is taken into account, incorporating the formation and destruction of hydrogen 
molecules, the radiative transfer is expected to be of greater significance (Haiman, Rees & Loeb 1996, 
1997). We will investigate the dynamical evolution of the gas clouds in this regime in future publications 
(Kitayama et al. in preparation; see also Susa & Umemura 1999). 

Our calculations are in good accordance with those of Thoul & Weinberg (1996), for the the same set 
of parameters, i.e. in the optically thin case and Zc = 2 — 5. This gives an important cross-check of the 
current results and confirms that they are insensitive to the different choices of initial conditions and central 
boundary conditions between their calculations and ours (see §2.4 and §4.1 of Thoul & Weinberg 1995). 

Based on the results of numerical simulations, we have predicted the global production rates of cold 
gas, hot gas, and unbound gas in the standard CDM universe. The abundance distribution of the cold gas 
matches well the observed mass ranges of dwarfs and galaxies. The global production rate of cooled gas is 
found to rise steeply from the local universe to ^ ~ 2 — 3, indicating a higher efficiency of star formation 
activity at high redshifts. 

We further predict that a large number of the Hell and Hel forest lines may arise in the quasar spectra 
at the UV wavelengths, which can be detected by HST and future space missions. Such observations should 
provide powerful probes of the physical state of the intergalactic medium, such as the gas kinematics, and 
the UV backgromid spectrum (e.g. Sargent et al. 1980; Ranch 1998). In addition, the existence of the 
warm neutral gas is inferred, due to the high efficiency of heating over ionization. The temperature of such 
gas can be as high as ~ lO^K and it may be related to the high spin temperatures suggested from the 21cm 
absorption line observations in the nearby damped Lya systems (Lane et al. 1998; Chengalur & Kanekar 
1999). In testing these predictions, it is essential to perform the multi-line analyses of the absorption 
systems, which will become possible with greater precisions in near future. 
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A. An analytical solution for radiative transfer: photoionization coefficients and heating 

rates of primordial gas 

We derive analytical formulae for the photoionization coefficients and heating rates of a radiation field 
penetrating into primordial gas aligned in a plane parallel geometry, with an arbitrary density profile. The 
frequency/direction-dependent radiative transfer due to the absorption by multiple species is explicitly taken 
into account. While the gas composed of atomic hydrogen and helium is considered here, the formalism can 
be readily extended to include other species. The derived formulae are applicable to a variety of problems, 
e.g., the reionization of the universe, and the photoionization of a mini-pancake or low-metallicity gas in 
galactic haloes. 



A.l. Derivation 

Suppose that an isotropic incident radiation field is processed through a gas slab composed mainly 
of atomic hydrogen and helium. For simplicity, we neglect the emission and scattering of photons by the 
intervening medium and only consider the absorption above the ionization energy of each species. We model 
the incident radiation spectrum as a power law with an index a: 

loiiy) = loimi) ( —) , (Al) 

where Iq^vui) is the specific intensity at the ionization frequency of neutral hydrogen j/hi. Then the 
processed mean intensity at an arbitrary point inside the slab is written as 

= ^J diohiv^i) (^-^^ (A2) 

using the optical depth t{i',uj) along a photon path s from the slab boundary to the point in the incident 
direction lj: 

t{v,uj) ^"^aiiv) I mds, (A3) 

where ai{i') and rii are the photoionization cross section and the number density of the species i (= 1,2,3 or 
HI, Hel, Hell, in the ascending order in its ionization energy), respectively. For simplicity, we approximate 
the cross sections by a single power-law: 

a^iiy) = aM) (^^) ©('^ " (A4) 

where Q{x) is the Heviside step function, and the amplitude ai{i'i), the index r/i, and the ionization energy 
hi'i are taken from Osterbrock (1989) and listed in Table |[ Note that the following formalism can be 
readily extended and applied as long as cross sections are expressed as a superposition of power-laws. 

Using the above expressions, the photoionization coefficients and heating rates for the species j are 
written as 

r MM 

= /..r^riiV°'"e-K-), (A5) 
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^3 = 



AnJiv) 



hv 

cJj{vj)Ig{vj) 



aj{i'){hv — hvj)dv 



duj 



dv f V 



-Vj 



(A6) 



Note that Tij defined above is related to the heating rate per unit volume used in the main text (eq. |^) 
hy TL — "^njUj. In general, t(v,uj) is a complicated function of v as it is the sum of components with 
different spectral indices (eqs |]A^ ] , [ A4 | ) . However, the frequency integrations in equations (M) and ( A6) 
are separated into intervals as 

3 



dv 



E 



dv with = oo. 



and it will not be a bad approximation to adopt for each interval 

-nf 



t{v,U}) ~ t{v,,uj) 



at Vj < V < V. 



(A7) 



(A8) 



where T{i'i,uj) takes the sum over all species at Vi (eq.|A3]), while the effective index 77°^ is set equal to the 
index 77 of the species which makes the dominant contribution to T{vi,Lo). Then the frequency integrations 
in equations ( [A5| ) and (A6) are performed analytically to give (see also Tajiri & Umemura 1998; Susa & 
Umemura 1999) 



7(/3ji,f(fi+i,a;)) / 



■'i+l 



(A9) 



with 



/3: 



nOff 



_ a + Tjj — 1 



(AlO) 



(All) 



where f^Vi+ijiu) is an approximated optical depth computed by equation (A8) at its upper limit Vi+i, and 
7(0, x) is the incomplete gamma function: 

-f{a,x)= [ e-H^-^dt, a>0. (A12) 
Jo 

In most astrophysically interesting cases, /3 > and /3' > are both satisfied. 

For a sufhciently long slab in a plane parallel symmetry with an arbitrary density profile, the optical 
depth given in equation (A3) varies with a direction cosine fi ~ cos0 (0 — if perpendicular to the plane) as 

t(i/,cj(i/)) = -^ai{v)N„ 



-r_i(i/) 



(A13) 
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where Ni = ^g_^nids is the column density perpendicular to the plane. Substituting equation (A13) into 
equations (|A^), ( A9), ( AlCl| ), and taking J dto 2tt dfi, we obtain the following analytical formulae for 
the photons propagating from one of the boundaries: 



J{u) 

r, 



- 



-f{Pji,T±{Vz+l)) 



Hj = hvj 



2TTaj{iyj)Ia{iyj) 



where f±(ui^i) = T^(^'i)(fi+i/!^i) , and a function f{a,x) is defined by 









H 







f{a,x) 



1 



a + 1 
( 1 



l{a,x) 



+ e - xEi [x) 



(a:^0), 

a+1 x'^ ^ ' 



using the incomplete gamma function (eq. | A12|) and the exponential integral: 



E^(x) = x" 



dt, a;>0, n = 0,1,2. 



(A14) 



(A15) 



(A16) 

(A17) 
(A18) 

(A19) 



The above results show that one can explicitly compute the photoionization coefficients and heating rates, 
once given the column densities perpendicular to the plane. In the optically thick limit t±_ — > oo, the 
photoionization coefficients and heating rates vary as (see also eqs [ A24 l ||A25t| below). 



A. 2. Simple example 



pure hydrogen gas 



If the gas is dominated by a single species, say atomic hydrogen, the expressions derived above can be 
reduced to even simpler forms as presented below. Such formulae are also quite useful in making physical 
estimations in many astrophysical problems. 

For the pure hydrogen gas, the effective indices rj^^ can be all set equal to 77^^ = ?7hi, and equations 
ij) and ( |A16| ) reduce to 

27ra-Hi(t'Hi)/o(j^Hi) /(/3hi,t_l(i'hi)) 



Thi 



mi 



hvm 



27rcrHl(i^Hl)/o(i^Hl) /(/3hi,t_i(i^hi)) 



mi 



HI 



(A20) 
(A21) 



with r^Hi — 3, /3hi — {a + 3)/3, and /3^j — (a + 2)/3. Given a specific value of a, these expressions are 
computed using equation ( |A18| ) in the optically thin and thick limits as 
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1. optically thin limit: Tj^(i^Hi) ^ 



HI 



Ti-m — 



2. optically thick limit: r^(i^Hi) 
Thi = 



oo 



1.49x10-12/21 s-i (a=l), 

7.47x10-13/21 s-i (a = 5), 

1.08x10-23/21 ergs-i (a = 1), 

2.33x10-24/21 ergs-i (a = 5), 



7.62 X 10-i3/2ir-^/3(t.Hi) s 



(A22) 
(A23) 



8.18 X 10-i3/2ir 



-8/3 



(a = l), 
(a = 5), 

2.17 X 10-23/2irji(i/Hi) erg {a = 1) 

1.55 X 10-23/21x7'^/^ 



(i^Hl) 

■hi) 

(f^Hl) 



erg s 



-1 (a = 5), 



(A24) 
(A25) 



where I21 = /o(t'Hi)/(10~2ierg s^^ cm-2 str-^ Hz-^). Note that equations (A20)- (A25) take account of 
the photons propagating from only one of the boundaries. In particular, the expressions in the optically 
thin limit (eqs A22 | |A23| ) should be multiplied by 2 to incorporate the photons coming from all directions. 



From the above results, one can readily estimate the ionizing structure of a medium exposed to the 
external photoionizing flux. We restrict our attention to the gas of temperature T ^ lO^K and ignore the 
collisional ionization. Then the ionization balance is expressed as 



Thi-'^hi ~ aH'^H(l - -'^Hl)^ 



where 



an ^ 3.96 x IQ-^^T^ 



"■^ cm^s-i, 



(A26) 



(A27) 



is the hydrogen recombination rate to the ground level (Spitzer 1978; Fukugita & Kawasaki 1994), and 
T4, = T/IO^K. Inserting equations ( |A22D and ( |A24D into (^2^), we obtain 



1. optically thin limit: t_l(z^hi) 

-'^HI 



(1 - ^Hl)2 



0.133 I^imT^" '^ (a = 1), 
0.265 I^imT^" '^ [a = 5), 



(A28) 



2. optically thick limit: T^(i^Hi) - 



(1 - Xm)' 



0.520 I^imT^^-^rl^^ii^m) 

0.484 /2->ir4-"-V^/^(i.Hi) 



(a = l), 
{a = 5), 



(A29) 



where ni = nn/cm and we have multiplied equation (A22) by a factor 2 to incorporate the photons 
coming from all directions. 
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Fig. 1. — Trajectories of radii of gas shells enclosing 0.2 (inner-most shell), 10, 30, 50, 70, 90 % of Mbound 
in the low-redshift collapse {zc — 0.5) for a cloud with Vc = 32 km s~^ (Mdoud — 9 x 10® M0) and different 
parameters of the UV background; (a) J21 = 1,^ = 1, (b) J21 = 0.1, a = 1, (c) J21 = 1, a = 5, and (d) 
evolving J21, a = 1. Different lines indicate the no UV case (dashed), the optically thin case (dotted), and 
the case with absorption (solid). 
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Fig. 2. — Evolution of temperature (highest panels), X^ii (second highest), TVhi (second lowest), and XhcII 
(lowest) of gas shells enclosing 0.2% (thick lines) and 90% (thin lines) of Mbound in the low-redshift collapse 
{zc = 0.5) for a cloud with Vc — 32 km s~^. (a) J21 = 1, a = 1, (b) J21 = 0.1, a = 1, (c) J21 = 1, a = 5, 
and (d) evolving J21, a = 1. Different lines indicate the optically thin case (dotted lines) and the case with 
absorption (solid lines). 
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Fig. 3. — Radial profiles at 2: = of (a) hydrogen density, (b) temperature, (c) HI column density from the 
boundary, (d) Xm, (e) Xb_c1, and (f) XhoII in the low-redshift collapse {zc = 0.5) for a cloud with Vc — 32 km 
s~^. Lines correspond to different parameters of the UV background as shown in the figure (unless indicated 
explicitly, absorption is taken into account). 
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Fig. 4. — Same as Fig. except for a cloud with (a) Vc — 50 km s ^ (Afcioud — 3 x 10^ M©), and (b) 
Vc = 100 km (Afcioud ^ 3 X 10^° Mq), in the case of J21 = 1 and a = 1. 




Fig. 5. — Same as Fig. except for tlie middle-redsliift collapse {zc — 3) with Vc = 32 km 

(Meloud==^2x 10»Mq). 
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6. — Same as Fig. g, except for the middle-redshift collapse {zc — 3) with Vc = 32 km 
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Fig. 7. — Same as Fig. |l|, except for the high-redshift collapse {zc ~ 10) with Vc = 32 km s 

(Meloud=i5x IO^Mq). 
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Fig. 8. — The fraction of gas mass collapsed after the onset of the UV background zyv = 20 as a function 
of circular velocity. Each horizontal row traces the time evolution at t = 0.5tc, tc, 2tc of runs with Zc = 0.5 
(top panels), 3 (middle), and 10 (bottom). Lines and symbols correspond to different parameters of the UV 
background as shown in the figure (unless indicated, absorption is considered). 
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Fig. 9. — The critical circular velocity V^^*" (thin lines) and Vj^^'* (thick lines) of the collapse (see text for 
definitions) with absorption (triangles, solid lines) and in the optically thin case (circles, dotted lines), (a) 
J21 = 1, a = 1, (b) J21 = 0.1, a = 1, (c) J21 = 1, a = 5, and (d) evolving J21, a = I. Larger symbols in 
panel (a) indicate V^^^ in the case of zuv = 50. 
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Fig. 10. — The baryon mass corresponding to T^'^^'^ as a function of collapse epoch Zc (thick lines). Different 
parameters of the UV background are assumed as shown in the figure (unless indicated, absorption is 
considered). Also plotted are the Jeans mass with T = 10** K (thin dot-short-dashed lines) and the masses 
corresponding to 1,2,3(7 density perturbations (thin dot-long-dashed lines) in the standard CDM universe 
with Slo = 1, ^ = 0.5, fib — 0.1, and erg = 0-6. 
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Fig. 11. — Density-temperature diagrams from our simulations and analytical estimates. Symbols indicate 
the virial density and temperature of simulated clouds in which the cooled gas mass reaches Mcioud before 
tc + tdyn (circles), before the present age of the universe to (triangles), and otherwise (crosses); (a) J21 = 1, 
a = 1, optically thin case, (b) J21 = 1, a = 1, with absorption, (c) J21 = 1, a = 5, optically thin case, and 
(d) J21 — I, a — 5, with absorption. Lines show the analytic relations icooi = idyn evaluated at Zc (solid), 
icooi = to — tc at Zc (dashed), and Tvir — 7cq at Zta, with the same UV parameters as the simulations except 
for assuming that the medium is optically thin in all cases. 
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Fig. 12. — Mass distribution of the baryon density that collapses per Hubble time at (a) — 10, (b) — 3, 
and (c) Zc = 0.5 in the standard CDM model with flo — 1, h = 0.5, fJf, = 0.1, and as = 0.6. Lines indicate 
the cases of different UV parameters or of no pressure as shown in the figure (unless indicated, absorption 
is taken into account). 
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Fig. 13. — (a) The baryon density which cools (thick hnes) or becomes bounded (thin hnes) per Hubble 
time under the UV background, normalized by that which collapses in the absence of gas pressure. Lines 
specify different UV parameters as shown in the figure (unless indicated, absorption is considered). The 
three regions divided by thick and thin lines correspond to "cold gas + star" , "hot gas" , and "unbound gas" , 
respectively, (b) Production rates of "cold gas + star" , "hot gas" , and "unbound gas" for the evolving J21 
and a = 1. For both panels, the standard CDM model is assumed. 
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Fig. 14. — The global production rates of "cold gas + star" per year predicted in the standard CDM model 
(unless indicated, absorption is considered). Also potted for reference are the observed cosmic star formation 
rates (Madau et al. 1996, 1998; Lilly et al. 1996; Connolly et al. 1997), without corrections for the dust 
extinction, compiled by Totani et al. (1997) and Totani (1999, private communication). 
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Fig. 15. — Line of sight column densities through the cloud Nj {j =HI, Hel, Hell) as a function of impact 
parameter p for Vc = 32 km s~^ (left panels) and Vc — 50 km s~^ (right). Thick lines are the clouds with 
Zc = 3 viewed at z = 3, while thin lines those with Zc = 0.5 viewed at z 0. Choices of the UV parameters 
(all with absorption) are as indicated in the figure. 
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Table 1: Characteristic redshifts of simulation runs 







^taO 


2cO 






low-redshift collapse 


40 


5.4 


3 


1.3 


0.5 


middle-redshift collapse 


110 


16 


10 


5.4 


3 


high-rcdshift collapse 


300 


45 


28 


16 


10 



Table 2: The critical impact parameter pcrt (see text for definition) and the ratio of its square for the evolving 
J21 and a = 1 







32 km s-i 




50 km s-i 


species 


Pert [kpc] 


ratio of ■pl^^ 


Pert [kpc] 


ratio of vl^^ 


HI 


10 


1 


10 


1 


z = 3 Hel 


0.8 


1/156 


1 


1/100 


HeH 


20 


4 


40 


16 


HI 


40 


1 


100 


1 


z = Hel 


10 


1/16 


10 


1/100 


Hell 


100 


6.25 


100 


1 



Table 3: Same as Table || except for J21 = 1 and a = 5. 







Vc = 


32 km s^i 


Vc = 


50 km s ^ 




species 


Pert [kpc] 


ratio of Pert 


Pert [kpc] 


ratio of p^rt 




HI 


20 


1 


30 


1 


z = 3 


Hel 


10 


0.25 


15 


0.25 




Hell 


100 


25 


150 


25 




HI 


10 


1 


10 


1 


z = 


Hel 


7 


0.49 


7 


0.49 




Hell 


200 


400 


400 


1600 
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Table 4: Indices of column density distributions, Nf cx p " and dJ^/dNf cx (Nf ) ^ 



species 






32 km s^i 




50 km s-i 


j 




11 


/3 


n 


/3 


HI (iVPj < IQi^cm- 


') 


1.5 


2.3 


1.3 


2.5 


HI (iVPj > IQi^cm- 


') 


3 


1.7 


24 


1.1 


Hel 




3.7 


1.5 


18 


1.1 


Hell 




2 


2 


2.3 


1.9 



Table 5: Relations among Nf,,, N?, , and N?, „ in cgs units 





A^gj < IQi^cm-^ 


iVgj > IQi^cm-^ 


J21 evolving, a = 1 


^^Hcl 


10"(^Hi/10'^)'-' 


/VP - 


^ lO-^iVPj 




/VP 


lOTVPj 


A/P . 

^^HcII 


■^1018 


J21 = 1, a = 5 


A/P ^ 

^^Hcl 


10i3(A^Pi/10i4)i-5 


AfP ^ 






A^P - 


lO^iVPj 


AfP ^ 

^^HcII 


■^10i« 



Table 6: Photoionization cross section parameters in equation ( A4) 



species 


ionization energy 


amplitude 


index 


i 


hvi [cV] 




rii 


HI(1) 


13.6 


6.30 X 10-18 


3.0 


Hel (2) 


24.6 


7.83 X 10-1* 


2.05 


Hell (3) 


54.4 


1.58 X 10-1* 


3.0 



